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Abstract 

We study the three-dimensional random Z(2) lattice gauge theory with Higgs field, which has the 
link Higgs coupling ciSUS and the plaquette gauge coupling C2UUUU. The randomness is introduced 
by replacing ci —> —ci for each link with the probability pi and C2 — C2 for each plaquette with the 
probability p2 . We calculate the phase diagram by a new kind of mean field theory that does not assume 
the replica symmetry and also by Monte Carlo simulations. For the case pi = P2(= p), the Monte Carlo 
simulations exhibit that (i) the region of the Higgs phase in the Coulomb-Higgs transition diminishes 
as p increases, and (ii) the first-order phase transition between the Higgs and the confinement phases 
disappear for p > pc — 0.01. We discuss the implications of the results to the quantum memory studied 
by Kitaev et al. and the Z(2) gauge neural network on a lattice. 
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1 Introduction 



Randomness and/or disorders are involved in quite many physical systems and play an important role 
in various fields of physics. In some electron systems, randomness appears due to impurities and affects 
electric conductivity. Sufficient randomness may cause localization of electrons [T]. 

For a quantum computer, the effects of environmental noises (including thermal fluctuations, disorders 
by impurities, etc.), which let the system to decohere, should be reduced to perform quantum computation 
as one designed. Kitaev'^] proposed a fault-tolerant quantum memory and quantum computations that are 
based on the Aharonov-Bohm effect of discrete Zi gauge symmtery. The memory is on a two dimensional 
torus and modeled as the Z(2) lattice gauge theory on a three dimensional (3D) lattice with the third 
dimension for the imaginary time. It is argued that the quantum memory works well if this model system is 
in the Coulomb phase] (an oredered phase with small fluctuations of gauge- field) instead of the confinement 
phase (a disordered phase with large fluctuations of gauge field). After that, many studies on Kitaev's 
model have appeared[3ji. Wang et al.[?j studied accuracy threshould of the model by using the Z-i random- 
plaquette gauge model(RPGM) on a three-dimensional lattice. In the Z(2) RPGM, the (inverse) gauge 
coupling for each plaquette takes the values ±/3 with random sign; the probability to take /3 is given by 
I — p and the probability of "wrong-sign" — /3 is p. The main interest is its phase structure, i.e., the 
critical concentration pdP) which disinguishes the confinement phase and the Coulomb phase. Wang et 
al.[4] calculated pdP = oo) ~ 0.029. Ohno et al. 5J calculated pdP) for finite /3 and showed Pc{f3) takes 
a maximum value pdP) — 0.032 ^ 0.033 at /? ~ 0.4 ^ 0.5. (We shall explain this pdP) in detail as the 
special case Ci = of the present model in Sect. 4 using Figdl) 

In this paper we consider another random gauge model related to the Z(2) RPGM, the 3D random 
Z(2) Higgs lattice gauge theory. In addition to the usual gauge coupling on each plaquette (i.e., the 
coupling of four gauge fields on the links around each plaquette) as the Z(2) RPGM, the energy of this 
model contains the Higgs coupling on each link (a pair of Higgs fields at the nearest-neighbor sites couples 
through the gauge field on the connecting link) . Both the coupling constants of Higgs and gauge couplings 
have wrong signs with certain probabilities. The pure case of the model is known to have three phases, 
Higgs, confinement, and the Coulomb phase[5]. We study its phase structure, in particular, how the phase 
boundaries shift as the randomness is increased. 

The reasons why we are interested in this model are as follows: First reason is related to the quantum 
memory. Are there any real materials that can be a candidate of Kitaev's model? Many strongly correlated 
electron systems can be described by an effective gauge field theory [7]. For example, the low-energy 
behavior of s = 1/2 antiferromagnetic Heisenberg spin model is described effectively by the U(l) lattice 
gauge theory coupled with CP^ spins. jS^ However, most of such low-energy effective gauge field theories 



^The ordered phase discussed in the Z{2) model of Ref. is the Coulomb phase instead of the Higgs phase as we shall 
clarify later. Because of the discreteness of Z(2) group, the mass of gauge boson becomes finite even in the Coulomb phase. 
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contains gauge coupling of matter fields in addition to the pure gauge term. When the gauge group is Z(2), 
the corresponding lattice gauge theory is the Z(2) Higgs lattice gauge theory, which is just the model we 
are going to investigate. In Kitacv's model, the Wilson loops are used as the index of quantum memory 
as well as the order parameter of confinement-deconfinement phase transition. However, for the gauge 
theory including matter couplings, it is well known that the Wilson loop always obeys the perimeter law 
and cannot be used as an order parameter. Another nonlocal order parameters are proposed [9 . We expect 
that this Higgs system is more realistic and can work as a quantum memory when it is realized in the 
Higgs phase. 

The second reason is related to neural network models of the human brain. [TOj The Hopfield model[TT] 
is well known as the standard neural network model to explain the mechanism of associative memory, i.e., 
a process of retrieving a pattern of neurons once stored in the brain. For the process of learning a pattern 
itself, various models based on the plasticity of synaptic connections have been proposed. In Refs.[6|,fT2]. we 
introduced and studied a neural network with gauge symmetry, in which the synaptic variables are treated 
as a gauge field. This sounds natural because electromagnetic signals propagate through the synaptic 
connections, and the electromagnetic interaction has U(l) gauge symmetry. There is a strong correlation 
between the performance of learning a pattern and retrieving it and the three phases of gauge dynamics. 

In neural networks, there should be necessarily malfunctions of signal propagations. In Ref. [12j. they are 
described as thermal fluctuations at finite pseudo temperatures by employing the Metropolis algorithm as 
the rule of time evolution as in the Boltzmann machine. Another possibility of description of malfunctions 
may be introducing impurities in the network. When the network is put on a lattice and the synaptic 
connections are restricted to the nearest-neighbor ones, such a model becomes just the random Z(2) Higgs 
lattice gauge theory we are going to study. It is interesting to see how the randomness caused by impurities 
affects the functions of human brain. 

The present paper is organized as follows. In Sec. 2, we introduce the model. In Sec. 3, we set up a mean 
field theory for systems with randomness and applies it to the present model. As the mean-field theory, 
the replica method for quenched averaged systems like spin glass is well known. The present mean-field 
theory is quite different from the replica method. It is general and applicable for any random system 
having "wrong-sign" coupling constants as long as usual mean field theory for the pure sustem is available. 
In Sec. 4, we study the phase diagram by Monte Carlo simulations. In Sec. 5, we preent discussion and 
conclusions. 

2 The Model 

The model is defined on the three-dimensional cubic lattice of the size V ~ L x L x L with the periodic 
boundary condition. On each site x there sits a Z(2) spin variable Sx — ±1 and on each link {x,x + /^) 
(/I = 1,2,3 is the direction index and we use it also as the unit vector) there sits a Z(2) gauge variable 
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Uxfi = ±1. The energy (multiphed by minus of the inverse effective temperature) A of the model is given 

by 

3 

A ~ ^ ^ ^ ^ Cixfj,Sx+fj,UxfiSx + ^ ^ ^ ^ C2x^j.vUxvUx+v.jj,Ux+^j..uUxfi, (2-1) 
an K /i<iy 

where ci^;^ and C2xiiv are random coefficients on each undirected link (x, a; + /i) and unoriented plaquette 
(a;, x + /i, x + /i + I', a; + i^, a::)(/x < v), respectively. The energy A of (|2.ip is invariant under the local gauge 
transformation, 

Sx^S'x^ WxSx. Uxf. ^ C/^^ = Wx+^.Ux^J.Wx. Wx = ±1. (2.2) 
cia;^ and C2a;^i/ are independent random variables taking the values as 



Cl with the probability 1 — pi 
— Cl with the probability pi 



C2 with the probability \ — P2 
C2x^u = { ^ _ (2.3) 

— C2 with the probability p2 

We regard the link with the "wrong sign" , ci^;^ = — ci, a link with impurity, and similarly, the plaquette with 
C2xtiv = —C2 a plaquette with impurity. Therefore pi and p2 are the concentrations of link and plaquette 
impurities, respectively, ci and C2 are positive parameters appearing in the pure system {pi = P2 = 0). 
We noite that, if one sets ci = then the present model reduces to the Z(2) random plaquette model 
considered in Refs.[H[l]. Each sample lattice has a fixed configuration of ci^,^ and C2xfiv To calculate a 
physical quantity like the internal energy E and the specific heat C which are measured in experiments, 
we first consider the thcrmal(anncaled) average {0{S, U)) of observable 0(5, U) for each sample as 



s u 



= ^^exp(AO =exp(-FO, (2.4) 

s u 

where i is the suffix to specify the sample as the «-th sample, and Zi is its partioton function and Fi is its 



free energy. Then we take the sample average {0{S, U)) of (0(5, U))i, a quenched average over sufficiently 
larage number Ns of samples, 

, Ns 

(0(5,C/))^ hm _^(0(5,J)},. (2.5) 

i—1 

This expression corresponds to physically measured quantities in a random system. The sample average 
(12.51) can be rewritten as 



(0(5,C/))= / [dci][dc2]pi2(c)(0(5,J))(c), (2.6) 
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where we wrote (0(5', J))(c) = (0(5", J)),,, [dci][(ic2] = Ylx,ij.'^^i^p-Ilx,ij.<,^ '^^^xt^'^^ introduced the 
probabiUty distribution of ci^fi and C2x^iv, 

Pl2{c)=W[{l-pi)5{cixii-Ci)+pi5{cixfi+Ci)\ J]^ [{I - P2)5{c2x,iu - C2) + P25{c2xfiu + C2)]- (2.7) 



The free energy the internal energy E per site, and the specific heat C per site are defined by 



3 Mean Field Theory for a Random System with Wrong-Sign 
Couphng Constants 

UsuaUy, the mean field theory (MFT) applied for a statistical system gives a rough but intuitive under- 
standing of the global properties of the system like its phase structure. We first summerize the MFT for a 
pure system and then develope a new MFT for a random system based on that for a pure system. 

3.1 MFT for a pure system 

For pure models without impurities a mean field theory can be formulated by the variational principle |13j. 
Let us briefly summerize it. We start with the partition function Z for the action A{(j)) with a set of 
variables 0, 



where F is the free energy. In MFT one prepares a trial action Ao(i^, A) having (a set of) variaitonal 
parameters A, the partition function of which is calculable; 



X,fL<l^ 



F = F(c), E 



{A){c)/V, C^{{A-{A)mc)/V. 



(2.8) 




(3.1) 




(3.2) 



Then there holds the following Jensen-Peiels inequality; 



F < F,{X) = Fo{X) + (AoiX) - A)o, 




(3.3) 



Thus, one minimizes Fy(X) by adjusting A to obtain the best approximation for F. 

For the present model ()2.1|) at pi = P2 = 0, it has been applied in Ref. 6 with the choice 




(3.4) 
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where ag, ai are the two variational parameters. Exphcitly, we have 

Fy/N ~ — log[2cosh(ao)] — dlog[2 cosh(ai)] + oq tanh(ao) + dai tanh(ai) 

— dci tanh^(ao) tanh(ai) ^ — - — tanh'*(ai). (3-5) 

The stationary conditions read 

flQ = 2dci tanh(ao) tanh(ai), 

ai — dci tanh^(ao) + 2(i((i — l)c2 tanh^(ai). (3-6) 

It predicts the three phases as hsted in Table In Fig[T]we plot the phase boundaries determined by 
MFT in the C2 — ci plane. 



phase 




{Sx)o 


ability 


Higgs 


^0 


^0 


learning and recalling 


Coulomb 


^0 





learning 


Confinement 








N.A. 



Tablel. Three phases predicted in MFT and the associated ability of neural net in a process of learning a 
pattern of Sx and retrieving it|12). 

3.2 MFT for a random system 

For the case of random system that involves quenched averages like ()2.5|) . one should generalize the MFT 
(|3.3p for a pure system. Below we develop such a generalization. For this purpose we start with the 
following variational system, which is independent of differences among samples; 

^o(a) = ^5]exp[Ao(5,i7,a)] =exp[-^^o(a)], 

aoxSx + Y^Y^aixtj.Uxfj,, (3.7) 

X X ^ 

where aox and aix^ are local variational parameters on the site x and link {x, x + ji) respectively. (We 
assign the suffices 0, 1, and 2 for site, link, and plaquette objects respectively.) By applying the inequality 
p.3p . one has 

F,{c) < Fy{a,c) = Fo{a) + (Ao)o(a) - (A)o(a,c). (3.8) 

^Some cases of Table 1 show the averages of gauge- variant quantities are nonvanishing, which violate Elitzur's theorem |14l . 
However, an additional averaging procedure over gauge rotated copies of the result of MFT gives rise to vanishing averages 
without modifying the phase structure of MFTfTS). 
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By multiplying pi2(c) both sides and integrating over cix^, C2xfii' one obtains the inequality for the free 
energy F, 



F< J [dci][dc2]pi2{c)F,{a,c). (3.9) 
Below we treat and ai^fj, as random variables described by the probability distribution poi(a) defined 



as 



i(a) = - qo)S{aox - oo) + Qo^ia-ox + Oo)] - qi)S{aixf, - ai) + qiS{aixf, + Oi)], (3.10) 



Po 

We consider the sample average of a function /(a) over samples of variational systems, each of which has 
different set of aox, o-ixfi, and write it as /(a), 

fia)= J[dao][dai]poi{a)fia). (3.11) 

By multiplying poi{a) both sides of (|3.9p and integrating over aox, oix/^, one has an upperbound Fy for the 
free energy F as 

F<Fy= j[dao][dai][dci][dc2]poi{a)pi2ic)F4a,c). (3.12) 

To relate pi,P2 and qo,qi we regard go,?! as concentrations of effective impurities on the sites and links 
respectively, and compose real impurities as the products of them. Explicitly, we determine the two 
functions 9o(pi,P2) and qi{pi,P2) so that the following relations hold in average over samples with poi(a); 

(1) average number of links for which signlaoxaixfiCLo.x+fj.] = — 1 is ^Vpi 

(2) average number of plaqucttes for which sign [a -1 is 3Vp2. (3.13) 

In Figl2]we present gi, 2(^1,^2) calculated by using a 3D lattice of the size 24^. 

Then Fy of ()3.12|) is calculable because J[dao][dai] and /[dco] [dci] are straightforward due to the 
decoupled nature of Fq, Aq and A. 

Fv/V = — log[2 cosh(ao)] — c?log[2 cosh(ai)] + ao tanh(ao) + dfli tanh(ai) 
-rf(l - 2po)(l - 2go)^(l - 2gi)ci tanh^(ao) tanh(ai) 



d{d~l 



^(1 - 2pi)(l - 2gi)''c2tanh''(ai). (3.14) 



2 

Fy becomes a function of two parameters ao,ai, which we adjust to minimize it. It is straightforward to 
see that the conditions of minimization take the same form as the conditions (|3.6p for the pure system 
{pi — P2 — 0); the former is given from the latter (the pure system) with the replacements 

Cl c[ = Gi{pi,p2)ci, 

C2 4 = 6*2(^1, P2)C2, 

Gi(pi,p2) = (l-2pi)[l-2go(pi,P2)]'[l-2gi(pi,p2)], 

G2{pi,P2) = (l-2p2)[l-2gi(pi,p2)]^ (3.15) 
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Thus the phase transition points (c^andom^ ^random ^ ^^le random system {pi^P2 0) are obatined by using 
the transitoin points (c^^^^^, C2^^^^) of the pure system {pi ^ P2 = 0) as 

1 mirr' 



„random 



(Pl,P2) = 



Gl(pi,P2) ^ 



1 



Because Gi^2(pi,P2) > 1, one obtains the following general inequalities 

random _ A \, pure random _ \ \ pure /f) ^ ^\ 

Ci (pi,P2)>q , C2 (pi,P2)>c5 . (3.17) 

This implies that the effect of impurities/disorder increases the critical coupling constants, which accords 
with our intuition. In Fig|3]we plot the critical lines of the present model in the C2 — ci plane obtained by 
(|3.16[) for the cases of p {= pi = P2) — 0.015 ^ 0.050 together with the critical lines for p — 0. 

We note that the present MFT can be applied for a random lattice system similar to the present one 
as long as its variational free energy is decoupled to single-site and single- link integrals like (|3.7p . In 
particluar, the probability distributions 9o.i(pi,P2) of effective impurities are universal, i.e., independent 
of each model, and determined by (|3.13[) . Therefore, for each critical point of the pure system, there may 
be one associated critial point of the random system as a modulation due to randomness like 



4 Monte Carlo Simulations 

In this section, we start with summerize the known results for pi = P2 ~ and for ci = 0. Then we 
study the phase struture of the 3D Z(2) random Higgs lattice gauge theory by MC simulations for the case 

Pl^P2= P- 

4.1 The case of p = 

Let us first summerize the phase structure for the pure case p = 0. The points in Fig[T] is the phase 
boundary in the C2 — ci plane [6], exhibiting three phases. 

The confinement-Higgs transition is first-order and terminates near (c2,Ci) ~ (0.55,0.35) as the com- 
plementarity argument predicts. 17J In fact, at C2 = 0, the partition function is exactly calculable by the 
single-link sum over Uxfi as 

Z|e,^o = [2cosh(ci)]3^, (4.1) 

which has no singularity in the ci-dependence and so there is no phase transition along C2 = 0. The 
Coulomb-Higgs transition is second order. In the limit C2 00, U^fj, becomes a pure-gauge configuration, 
Uxfj. = Wx+^Wx {Wx = ±1), the system reduces to the 3D Ising spin model A — ci J2x ^ ^'x+^^'x i^'x = 
WxSx), which exhibit a second-order phase transition at ci ~ 0.22. The confinement- Coulomb transition 



8 



is also second order. Along ci — 0, the system reduces to the pure gauge model, which is known to exhibit 
a second-order phase transition at ci ~ 0.76. 

4.2 The case of ci = 

Next, let us summerize the case of no Higgs coupling (c\ — 0) but with randomness. In Fig|4]we present 
the phase diagram in the p — T{p = p2, T = C2^) plane obtained in Ref.[S]. At p = 0, a second-order phase 
transition at the critical point T = T^. ^ 1/0.76. For T < the gauge-field fluctuations are small and the 
system is in the ordered Coulomb phase, while for Tc < T, the fluctuations are large and the system is in 
the disordered confinement phase. As p increases from p = 0, the critical coupling constant Tc{p) decreases 
due to the randomness in gauge couplings, and it shall vanishes at a certain value po, Tc{po) = 0. The 
data of specific heat suggests that the order of transition changes from second order to higher order for T 
smaller than about 1.15. Tdp) seems to show a "reentrance behavior", i.e., as T is lowered, the value of p 
on the curve Tc{p) increases first and takes the maximum value Pmax — 0.033[5j around T ~ 0.4 ~ 0.5 and 
then decreases to end up with p ~ 0.029 at T = 0[4]. 

4.3 Set up of MC simulations and the global phase structure 

In our MC simulations we use the standard Metropolis algorithm [TS]. The typical number of sweeps in a 
single run for {0)i of each sample is ~ 2 x 10^. To estimate the error of {0)i, which we call the MC error, 
we divide a run into 10 successive intervals to generate 10 data. AS the MC error we estimate the standard 
deviation of these 10 data. The average acceptance ratios are 0.4 ~ 0.5. In FigEK we present a typical 
result of the specific heat C together with the MC errors. For a random system one may be interested in 
the standarde deviation over samples (SDS), which should converge to a nonvanishing value even in the 
limit Ns oo. In FiglSb we present such SDS for the same C as FiglSK with the sample number Ns = 40. 
By comparing these two figures, both quantities are similar in magnitude (they differ by up to factor ^ 2), 
and have similar behaviors (i.e., the ci dependence). In the figures below the error bars show SDS. 

To see the dependence of (O) on Ng and to find a suitable value of Ng we present in FiglHthe specific 
heat C vs Ng. It shows that the Ns dependence is rather weak for Ns > 20. Below we show the results 
for Ns — 40 otherwise stated, where the error bars in the figures show SDS. 

Let us first observe how the peak of the specific heat shifts as p increases. In FigllK we present the 
curves of the peak location of specific heat in the C2 — ci plane. In FiglTb, we present the specific heat for 
various p at C2 = 1.0. As p increases, its peak becoms rounder and its position shifts to larger ci direction. 
It is consistent with the MFT prediction of Figl3] in Sect. 3. In Fig[7l;, we present the specific heat for 
various p at ci = 0.1. Its peak becomes rounder but its position is almost unchanged. We shall discuss the 
new behavior of C appearing in larger latice systems. 
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4.4 Study of the three cases 

Let us see whether these specific-heat peaks exhibit genuine phase transitions or not. To be specific, we 
focus on the following three cases: 

(A) The confinement-Higgs transition along C2 = 0.7, 

(B) The Coulomb-Higgs transition along C2 = 1.0, 

(C) The confinement- Coulomb transition along ci — 0.1. 

(A) C2 = 0.7: 

For p — Q we know that this case gives rise to a first-order phase transition. In FiglSl we present U and 
C for p — 0.01. U shows no hysteresis and the peak of C shows no systematic system-size dependence, 
which imply that this peak implies a higher-order transition or a crossover but not a first-order nor second- 
order phase transition. The small rendomness of p = 0.01 is sufficient to destroy the first-order transtion 
of p = here. Because the firsr-order transition line for p = ends at C2 ~ 0.55 and no phase transitions 
follow in the smaller C2 region (see eq. (|4.ip ') it is quite possible that there are no genuine phase transitions 
of finite order in the case (A) for p > 0.01. 

(B) C2 - 1.0: 

For p = this case gives rise to a second-order phase transition (See FigH]). For p = 0.01, 0.02, 0.03, 0.04 
the specific hear shows systematic size-dependent development, which supports that the second-order 
transition survives up to p = 0.04. The peak of C for p ~ 0.05 fails to show systematic size dependent 
development. In Figl9^,b, we compare C for p = 0.04 and p — 0.05. Even for the size up to L = 16, the 
difference appears clearly. To confirm the existence of second-order transition for p = 0.04, we fit C of 
p = 0.04 with L — 24, 28, 32 shown in Figl^^ by the finite-size scaling hypothesis ^20^ which reads for C as 

C(cuL) = L^ f{Lie), , ^ fl^fl^, (4.2) 

where f{x) is the scaling function, cioo is the critical value of ci for the infinite system L — oo, and a, v 
are the critical exponents. In FiglHb we present /(x) calculated with the choice 

cioo = 0.277, V = 0.90, a = 0.273, (4.3) 

which supports the existence of /(x). For p = 0.05 such a fit was impossible due to the lack of systematic 
size dependence. 

(C) ci = 0.1: 

For p = we know that this case gives rise to a second-order phase transition (See FiglT]). We shall 
argue that this transition changes to a third-order one for p — 0.01 and becomes a higher-order one or a 
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cross over for p = 0.02. In FigfTUb we present a close up of C near C2 = 0.80 for p = O.Oljf] It shows that a 
small peak is developed, which shifts, as L increases, to larger C2 direction and diminishes gradually. We 
interprete this behavior of C implies that the system is just in the transient point from the second-order 
transition to a higher-order one or to a crossover. We remember that the possibility of similar change from 
a second-order transition to a higher-order one has been pointed out in Ref. 5 for ci = as one lowers 
the critical value of c^^ (See Figd]). To study this point in detail, we measured a new observable, the 
derivative dC/dT of C. To introduce the "temperature" T we multiply the action by a factor (3 = T^^ 
and making the derivative of C w.r.t. T, 

dC 



Sample average of 
dT 



P^iA") - {Af)+p\{A^) ~ (A)2)(3/3(A) -2) /V, (4.4) 



and set /3 = 1 finally. In FigllOb we present dC/dT. Its signature changes at C2 ~ 0.83{L = 32) where the 
small peak locates. This may be a precursor of a third-order transition. It is natural that this small peak 
disappears at L — > oo and a sharp edge (cliff) remains, which implies a thrid-order transition characterized 
by a finite gap A{dC/dT) ^ 0. In Figs[lOl:,d we present C and dC/dT for p = 0.02. There are no 
indications of possible third-order transition. They may describe a smoother transition or a crossover. The 
detailed study is necessary to calculate pc of the present case, which may involve third-order or higher-order 
phase transition. As explained above, this is consistent with the result of Ref. [5]. 



5 Conclusion and Discussion 

In this paper, we considered the 3D random Z(2) Higgs lattice gauge theory and studied its global phase 
structure in the C2 — ci plane by a new MFT and MC simulations. The MC simulations showed that, 
for the case (A) of confinement-Higgs transition, the first-order transition for p = disappears quickly at 
p = 0.01. For the case (B) of Coulomb- Higgs transition, the second-order transition for p = persists up to 
p = 0.04. The result for p — 0.05 is failed to fit the scaling law. For the case (C) of confinement- Coulomb 
transition, as L increases, a small peak in C developes for p = 0.01 and then diminishes, which leads to 
a third-order transition. However, for p = 0.02, a third-order transition has not been observed, so further 
study is necessary to explore a possible higher-order transition. 

Let us discuss some implications of these results for a quantum memory. Inclusion of the Higgs coupling 
produces the third Higgs phase which has the least fluctuations of variables among the three phases, and 
so more stable functions than in the Coulomb phase are expected. The estimated critical concentration 
Pc — 0.04 ~ 0.05 for the case (B) is larger than pc ~ 0.033 for ci = 0. Of course, the nonlocal observables 
become more complicated than the Wilson loop itself [9]. 

Concerning to the neural network, "thermal" fluctuations previously considered as noise effects in 

^By comparing Fig llOb with C for ci = 0[5] we find the two peaks of C for L = 24 locate at almost the same place 
C2 ~ 0.82. 
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Ref . [51 [T2] act uniformly in space. In contrast, the random effects considered in tlie present model are 
inhomogeneous. However, these two effects seem to have no crucial qualitative differences for a system 
staying near a phase boundary but just inside the Higgs phase. Either introduction of tiny amount ofran- 
domnessin the local coupling constants or heating the sytem by tiny amount of temperature rise necessarily 
drives the system into the Coulomb or confinement phase. 

When one tries to include the time evolution to a 3D neural net in a manner faithful to quantum theory 
and/or statistical mechanics, one faces a 4D system (with the imaginary-time as the fourth direction). 
Then the effect of randomness to such a 4D system in of interest. Some analyses have appeared for 4D 
Z(2) RPGM[5T]. Including the Higgs coupling to this 4D model may be an interesting subject to study. 
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Figure 1: Phase diagram by MFT for p = 0. The dots represents the results of MC simulations .[6] The 
dashed curve C2 — —^£nth{ci) is the exact results obtianed by the self-duality argument jl6|. on which 
some of the transition points may lie. 




Figure 2: qo{pi,p2) and qi{pi,p2) satisfying (|3.13[) calculated by using a 3D lattice of the size 2A^. (a) 
P2 vs qi for pi ~ 0.01,0.03,0.05. qi is independent of pi because the wrong p2-plaquettes are generated 
by gi-links alone without go-sites, (b) p2 vs qo for pi — 0.01,0.03,0.05 from below, qo decreases as p2 
increases because too sufficient number of gi-links are generated to make p2-plaquettes. (c) pi vs qi for 
P2 — 0.01,0.03,0.05 from below, qi is independent of pi because qi is fixed only by the number of p2- 
plaquettes. (d) pi vs go for P2 — 0.01,0.03,0.05 from above. For very small pi, gg-sites are unnecessary 
because of sufficient qi -links supplied to compose p2-plaquettes. 
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Figure 3: The phase diagram in the C2 — ci plane determined by the mean-field theory for p = 
0, 0.015, 0.022, 0.035, 0.05 from below. As the disorder increases, the region of Higgs phase diminishes. 
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Figure 4: Phase diagram in the p — T(p = p2,T — C2^) plane for ci = from Ref. fS^ . The phase boundary 
p = pc{T) gives the maximum value p ~ 0.0332 ~ 0.033 at T ~ 0.3 - 0.4. It ends at T with 
p ~ 0.029 4J. The solid curve is the Nishimori line [19] on which the "thermal" effects valances with the 
"random" effects. 
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Figure 5: Specific heat for C2 = 1.0 and p = 0.04 with L = 24, = 40. (a) With MC errors and (b) With 
standard deviation over samples (SDS). 
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Figure 6: Sample number {Ng) dependence of the specific lieat for C2 = 1.0, ci = 0.27 with p — 0.03 and 
L = 12 with SDS. 
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Figure 7: (a) Phase diagram in the C2 — ci plane with the possible phase boundary curves determined by 
the location of the peak of specific heat C with L = 12. The three segments marked by (A,B,C) show the 
cases we shall examine in detail below, (b) p-dependence of the specific heat C for C2 = 1.0 with L = 12. 
(c) p-dependence of the specific heat C for ci = 0.1 with L = 12. 
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Figure 8: Inernal energy U and the specific heat C for the case (A) C2 = 0.7 at p = 0.01. (a) U has no 
hysteresis and (b) C has no systematic size dependence. 




Figure 9: Specific heat C for the case (B) C2 = 1.0. (a) C for p = 0.04; (b) C for p = 0.05. (c) C for 
p — 0.04 before scaling; (d) C for p ~ 0.04 after scaling. 
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Figure 10: C and dC/dT for the case (C) ci = 0.1. (a) C for p = 0.01, (b) dC/dT for p = 0.01, (c) C for 
p — 0.02, (d) dC/dT for p = 0.02. In each figure the SDS are shown separately. 
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